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We study a two-dimensional atomic mixture of bosons and fermions cooled into their quantum 
degenerate states and subject to an optical lattice. The optical lattice provides van Hove singu- 
larities in the fermionic density of states. We find that these van Hove singularities produce new 
and interesting features for the transition towards phase separation: an arbitrary weak interaction 
between the bosons and the fermions is sufficient to drive the phase separation at low temperatures. 
The phase separated state turns stable for attractive and repulsive interaction between the bosons 
and fermions and can be cast into the standard form of a 'liquid-gas' transition. 



I. INTRODUCTION 

Atomic mixtures of bosons and fermions represent a 
new laboratory system for the study of quantum degen- 
erate matter. A major experimental breakthrough is the 
sympathetic cooling of atomic fermions into their quan- 
tum degenerate state |^ |^ Q allowing for the realiza- 
tion of new thermodynamic phase transitions. The main 
focus is on the superfluid transition of atomic fermions 
0, 0, 0, and the appearance of complex quantum 
phases in the presence of an optical lattice, such as^super- 
solids 3 and nontrivial Mott insulating phases 10]. In 
addition, a Bose-Fermi mixture cooled into its quantum 
degenerate state shows a strong tendency towards phase 
separation and the demixing of the bosons and fermions 
12J. While this phase separation itself exhibits many 
fascinating phenomena [Tj, [Ij, [Tj, , it also represents 
a major drawback of Bose-Fermi mixtures as it restricts 
the allowed parameter range for the observation of com- 
plex quantum phases, see e.g., the competition between 
phase separation and a supersolid phase as discussed in 
Ref.H The analysis of phase separation in Bose-Fermi 
mixtures then contributes to our basic understanding of 
atomic gases as needed in future studies of other novel 
quantum phases. 

The main mechanism driving the phase separation 
is a small perturbation of the bosonic density 5n^ in- 
ducing a modulation of the fermionic density Sn^ ^ 
— UT,sN{eT,)5n^. Here, iV(eF) denotes the density of 
states at the Fermi energy and C/fb is the strength 
of the interaction between the bosons and the fermions. 
This fermionic distortion acts back onto the bosons and 
produces a shift of the bosonic energy —U^gN{e^)5n'^/2, 
thereby inducing an attraction between the bosons with 
strength U^gN{ep). Phase separation appears if the in- 
duced attraction is of the same order as the intrinsic re- 
pulsion Ub between the bosons. 

Previous work on phase separation in Bose-Fermi mix- 
tures concentrated on the T = limit. The phase sepa- 
ration is triggered by changing the interaction U^b and 
results in minimal energy configurations involving sepa- 
rated states with two pure phases, two differently mixed 
phases, or one pure phase coexisting with one mixed 
phase. Assuming a repulsive interaction between the 



fermions and the bosons, the phase separation of Bose- 
Fermi mixtures in a 3D harmonic trap results in an in- 
creased bosonic density accumulated in the center of the 
trap and surrounded by a fermionic shell 0, [T3 . [T3| . 
A thermodynamic analysis of the homogeneous situation 
has shown that the criterion for phase separation 
described above is preempted by a first-order scenario 
resulting in a separation involving a mixed and a pure 
phase. On the other hand, an attractive interaction be- 
tween the bosons and fermions t rigg ers a system collapse 
above a critical boson number [l6|; such a collapse has 
been recently observed in an experiment by Modugno et 
al. 

Here, we concentrate on the behavior of phase sep- 
aration in a situation, where the fermionic density of 
states is strongly modified by the presence of an opti- 
cal lattice. Of special interest is the two-dimensional 
(2D) situation where the presence of saddle points in 
the fermionic dispersion relation eF(k) provides loga- 
rithmic van Hove singularities in the density of states 
N{e) ~ iVoln|2WF/(e - eo)|, with the band width 
and eg the position of the van Hove singularity. A fi- 
nite temperature T cuts this singularity and the system 
response becomes strongly temperature dependent, re- 
sulting in a smooth transition with decreasing tempera- 
ture into a phase-separated state with two mixed phases. 
Below, we present a detailed analysis of this finite tem- 
perature phase transition. We find that the van Hove sin- 
gularities strongly enhance the tendency towards phase 
separation, with an arbitrary weak interaction between 
the bosons and the fermions already sufficient to drive 
the transition at low temperatures. The transition is of 
the 'liquid-gas' transition type and the phase separated 
state is stable for attractive as well as for repulsive inter- 
action between the bosons and fermions. Concentrating 
on the weak coupling limit, we derive the nT^-T phase dia- 
gram involving two mixed Bose-Fermi gases with slightly 
different fermion density Suy- Here, we focus on phase 
separation and ignore other possible instabilities in the 
system, e.g., BCS-superconductivity or supersolid forma- 
tion. Note that these competing phases also exhibit a 
strong increase of the critical temperature due to the 
van Hove singularity, a phenomenon well studied in BCS- 
superconductivity |l7| . Furthermore, our weak coupling 
analysis with U'^^Nq ^ U-b excludes the appearance of 
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a competing first-order phase transition of the type dis- 
cussed in Ref. jl^ 

The Hamiltonian of an interacting Bose- Fermi mixture 
subject to an optical lattice is presented in Sec. |n| We 
make use of the tight-binding approximation to cast the 
Hamiltonian in a form describing a Bose-Fermi mixture 
on a lattice. We then integrate out the fcrmions and ar- 
rive at an effective boson Hamiltonian. The interaction 
becomes strongly temperature dependent and the effec- 
tive boson Hamiltonian exhibits an instability towards 
phase separation at a critical temperature Tps. This in- 
stability is further examined in Sec. IHII and the phase 
separated state is studied within the Thomas-Fermi de- 
scription; the results are summarized in a phase diagram 
exhibiting all the properties of a standard liquid-gas tran- 
sition. Finally, we study the impact of this transition on 
an atomic Bose-Fermi mixture in a finite trap in Sec. II VI 



II. BOSE-FERMI MIXTURES 

We start with interacting bosons and fermions in two 
dimensions. Such a 2D setup is achieved via applica- 
tion of a strong confining potential in the transverse 
direction. The Hamiltonian for interacting bosons and 
fermions subject to an optical lattice takes the form 
H = Hs+Hp+ i?i„t with 
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-A + yB(x) Vb, 



(1) 



dx ( ^fbI/'b ^AbV'^V'f + -5bV'b V'b V'bV'b 



Here, ip^ and V'b denote the fermionic and bosonic 
field operators, while the interaction between the par- 
ticles is taken into account within the pseudopoten- 
tial approximation. We assume a repulsive interaction 
5b = 47raB?i^/mB between the bosons with the scatter- 
ing length Ab > 0, while the coupling g^^ = 27raFB^^/A' 
accounts for the interaction between the fermions and 
the bosons with /i the relative mass and Ofb the boson- 
fermion scattering length of either sign. Furthermore, we 
restrict the analysis to spinless fermions; such a spinless 
fermionic atomic gas is naturally achieved in an exper- 
iment via spin polarization. Then, the s-wave scatter- 
ing length of the fermion-fermion interaction vanishes, 
while p-wave scattering is suppressed at low energies 
and is neglected in the following analysis. The opti- 
cal lattice with wave length A provides an a = A/2- 
periodic potential for the bosons {a — B) and fermions 
(a = F) with V„(x) = Va sm'^{'Kx/a)+Va sm^{TTy/a)+V^, 
while the trapping potential accounts for the strong 
transverse confining rUaUj'^^z'^ /2 establishing a two di- 
mensional setup, and for the weak longitudinal trapping 



A. Hamiltonian within tight-binding 
approximation 

In the following, we focus on strong optical lattices 
ma and weak interactions. Then, 
in analogy to the mapping onto the Bose-Hubbard model 
T^, we transform the Hamiltonian ^ for the Bose- 
Fermi mixture to a simplified Hamiltonian of the tight- 
binding form. We restrict the analysis to the lowest 
Bloch band and introduce the Bloch wave functions Vk 
(fermions) and wj^ (bosons) of the corresponding single 
particle problem in the 2D periodic potential. In turn, 
the Bloch wave functions Wk(x) and tik(x) are related to 
the Wannier functions t5(x — R) and w(x — R) according 
to 



zz;(x-R) = ^ ^ w;k(x) exp (-zRk) , (2) 

kG-fsr 

w(x-R) = ^ i;k(x)exp(-iRk), (3) 
keif 

with R a lattice vector and K the first Brillouin zone 
of the reciprocal lattice. Here, we have introduced the 
quantization volume V — No? with N the number of unit 
cells. In the following, rip.B denote the number of parti- 
cles per unit cell. We express the bosonic and fermionic 
field operators V'f.b in terms of the Bloch wave functions 
Wk and Wki or equivalently, in terms of the Wannier func- 
tions w and V 



W - ^E^k(x)fe+=^zI5(x-R)fe+, (4) 

^ keK R 

P+(x) = ^E^k(x)c+=E^T(x-R)c+. (5) 



The bosonic creation operator of the Bloch wave state 
Wk is the Fourier transform of the bosonic creation oper- 
ator 6+ for the Wannier state {(;(x — R) with 



(6) 



and analogously for the fermionic creation operators 
and c+. In our subsequent analysis it is convenient to 
use the operators and . Inserting the expansion 
into the Hamiltonian and restricting the analysis to 
on-site interactions, the Hamiltonian of the Bose-Fermi 
mixture 1^ reduces to 

H = Y.e^{k)b+b^ + ^ E Wk'V (7) 



keif 



{k,k',q,q'} 



q£if 



N 



{k.k',q,q'} 



The summation {k, k', q, q'} is restricted to k, k', q, q' G 
K and the momentum conservation k — k' -|-q — q' = K™ 



3 



with Km a reciprocal lattice vector; a scattering event 
involving such a vector is an Umklapp process. The 
interaction parameters involve the Wannier functions w 
and V and take the form 

(7b = gnj dx|iJ(x)|4, (8) 
- / d^\w{^)\^\v{^)\\ (9) 



while (k) and (k) denote the lowest energy bands for 
the fermions and bosons, respectively. 

In a two dimensional system, the band structure eF(k) 
exhibits saddle points at wave vectors koi and energies 
eoi = eF(koi)- Lattice symmetries naturally produce sev- 
eral saddle points at the same energy eq; we denote their 
number by z. Close to such a saddle point, the dispersion 
relation eF(k) is quadratic, and can be expanded along 
the principal axes (k = kg + k_L + kn ) 



A reliable estimate of the interaction strength is pro- 
vided by approximating the Wannier functions z«(x) and 
^(x) by the wave function of the harmonic oscillator in 
each well. The oscillator frequency in each well is given 
by w^cii ~ y/^EaVa/h with a — B, F, which implies a size 
a^^i, = y^Ti/muj^^ii for the localized wave function. The 
interaction strengths [/fb and Ub, cf. © and ||5J), assume 
the form 



Ub 



C/f 



sV^Eb^ 

A 



VEbE^ 





Ti 


'Vb\ 


2Eb 




.Eb) 


Ofb / 




hLL>B± 


A I 


2E^ 


2Eb 



1/2 



1/4 



Vb 
Eb 



(15) 

1/4 



The validity of the derivation of the Hamiltonian (TJ re- 
quires that the interaction parameters Ub and [/fb are 
small compared to the energy gap ~ ?iw„eii separating 
the lowest Bloch band from the next higher. 



eF(k) eo 



2m I 



2m\i 



■II- 



(10) 



The saddle point gives rise to van Hove singularities in 
the density of states (per site) for e ^ eo, 



N (e) - Nq In 



2PFf 



e - eo 



(11) 



with the band 



width 

2\ 



z(mj^m||)i/2(a2/27r2?i^). 



Wp, and the prefactor Nq = 
This logarithmic singularity 
plays a crucial role in the study of phase separation in 
Bose-Fermi mixtures and will be examined in detail in 
the following. Note, that the situation is different in 3D 
systems, where the density of states is regular, and sin- 
gularities appear only in quantities involving derivatives 
of the density of states. 

As an example we study the dispersion relation for a 
strong optical lattice with only nearest neighbor hopping. 



eF(q) = -2Jf 



(12) 



and Jf the hopping energy. The density of states is ir- 
regular and exhibits a logarithmic van Hove singularity 
at eo = arising from saddle points at ko — (tt/o, 0) and 
ko = (0, 7r/a), 



N{e) = NaK 



'1- 



16J2 





16Jf 


- TVoln 




e 



(13) 



with iVo — l/{2n'^Jp), and K[m] the complete elliptic 
integral of the first kind . The hopping amplitude Jp 
derives from the exactly known width of the lowest band 
in the ID Mathieu equation 20] 



4Jf = 



16 



E^ 



1/4 



exp 




(14) 



B. Effective boson Hamiltonian 

In order to study the stability of the ground state, 
it is convenient to derive an effective Hamiltonian for 
the bosons alone. This effective Hamiltonian accounts 
for the fermions via a modified interaction between the 
bosons. We start with linear response theory, where 
the boson density nB(q) drives the fermionic system 
(?iF(q)) = UpBXiT,q)nB{q) with x(T, q) the response 
function of the fermions at temperature T. This per- 
turbed fermionic density in turn acts as a drive for the 
bosons and is accounted for by the effective interaction 
between the bosons 



1 



2N ^ 

{k,k',q,q'} 



[C/B + C/Flx(r,q-q')] hlh^,hlh^ 



(16) 

The response function of the fermions is given by the 
Lindhard function 



x(r,q) 



dk /[eF(k)]-/[eF(k-fq)] 
K vq eF(k) - eF(k + q) + z?7 



(17) 



with vq = (27r/a)'^ the volume of the first Brillouin zone. 
The temperature T enters via the Fermi distribution 
function /(e) = 1/{1 + exp[(e — ^p)/?"]}- We focus on 
static instabilities of the ground state and neglect the 
frequency dependence of the response function. A rigor- 
ous derivation of the effective action including its time 
dependence can be achieved within a path integral ap- 
proach ,21j . 



III. PHASE SEPARATION 

A. Instability at g = 

The Lindhard function is always negative and induces 
an attraction between the bosons independent of the at- 
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tractive/repulsive nature of the original coupling U-pu be- 
tween the bosons and fermions. The effective long dis- 
tance scattering parameter for q ^ takes the form 



(18) 



For a fermionic system with a regular density of states, 
the Lindhard function at q = and low temperatures re- 
duces to x(T ^ 0,0) = —N{ep) with N{e) the fermionic 
density of states and ep the Fermi energy. For fermions 
on a square lattice in 2D the density of states is irregular 
and exhibits a logarithmic van Hove singularity at eo, see 
Eq. (|13l) . For a fermionic filling such that the Fermi en- 
ergy e-p matches the energy of the van Hove singularity, 
i.e. e-p = eo, the Lindhard function diverges logarithmi- 
cally for T ^ 0, and its asymptotic behavior takes the 
form 



X(r^0,0)- / deN{e) 



. dm 

de 



-No\n 



2ciWp 
T 



(19) 



with ci = 2exp(C)/7r « 1.13 a numerical prefactor and 
C w 0.577 the Euler constant. 

As a consequence of this logarithmic divergence of the 
Lindhard function for T — > 0, the effective scattering pa- 
rameter U^it always turns negative at low temperatures. 
Since a thermodynamically stable superfluid condensate 
requires a positive effective interaction U^tt > the 
system exhibits an instability at the critical temperature 
Tps defined via ;7,ff(rps) = 0. Using Eqs. (O and (O, 
we find the critical temperature 



Tps = 2ciW^pexp [-1/Apb] 



(20) 



with ApB — {Up^/U^)Nq the dimensionless coupling con- 
stant. Note, that weak coupling requires Afb < 1 and 
the critical temperature Tpg is well below the transition 
temperature Tkt of the superfluid condensate. Below 
the critical temperature Tpg the effective interaction [/^f 
turns negative providing a negative compressibility, and 
the homogeneous superfluid condensate becomes unsta- 
ble. Then, the new ground state with fixed averaged 
density and Up exhibits phase separation with areas 
of increased and decreased local densities coexisting. 

Note, that this transition towards phase separation ex- 
hibits two major differences as compared to the phase 
separation discussed in Refs. First, the phase 

separation is an instability appearing at low tempera- 
tures for arbitrary small coupling Up^ between the bosons 
and fermions. Second, the increase/decrease in the 
bosonic density drives the fermionic density away from 
a filling factor close to the van Hove singularity provid- 
ing a regular x(2^j 0) which in turn stabilizes the system. 



B. Thomas-Fermi approximation 

In the following, we study the phase separated state 
within the Thomas-Fermi theory. Within this theory, 
we introduce two densities nB(x) and np{yi) which are 



smooth on a scale large compared to the Fermi wave 
length 1/fcp and the bosonic coherence length ^ — 
l/(87maB)^/^. Then, the syste m is in thermodynamic 
equilibrium at every position and neglecting the kinetic 
energy of the bosons the free energy T [tib , ^f] of the 
system becomes 

T[n^,np] ^ y"^|FpK(x)] -f V;*(x)np(x) (21) 



-|-Vb*(x)"b(x) + -?7BnB(x)nB(x) -|- Upsns{-x.)np{y.) 



Here, we include the weak external trapping potential 
Vb* for the bosons and Vp for the fermions. Furthermore, 
the appearance of the instability requires that the criti- 
cal temperature Tpg is large compared to the mean level 
spacing introduced by the external trapping potentials, 
i.e., Tps ^ hoj with uj a characteristic trapping frequency. 
Note, that we assume T <C Tkt which implies that the 
infiuence of thermally excited bosons is small and can be 
neglected. The local free energy Fp[np] of the fermions 
at temperature T takes the form (f = — ^n) 



Fp[np] = -T j deN{e) In 



1 + exp 



T 



/zrip 



(22) 

The local chemical potential /i(x) is determined by the 
condition that the fermionic free energy Fp is a minimum, 
i.e., d^Fp — 0, and hence 



(x) = / deN{e)f [e ~ ^^(x)] 



(23) 



Furthermore, the total number of particles in the system 
is fixed providing the additional constraints 



Np 



d'x.np{x), N-B 



dx TT-B (x) . 



(24) 



The ground state configuration at temperature T is de- 
termined by the minima of the functional J^[nj^,np] sat- 
isfying the conditions 12411 . and hence 

^K(x)]+V;*(x) + C/pBnB(x) = fip, (25) 
t/BnB(x) + VB*(x)-hJ7pBnp(x) = /Xb, (26) 

with Hp and the Laplace multipliers introduced to 
account for H24|l. These parameters act as global chemi- 
cal potentials for the system (we have used the relation 
d„pFp[np] = n[np]). 

In the following, we neglect the trapping potentials 
Vp and Vb* ^^'^ analyze the stability of the state with a 
homogeneous fermionic and bosonic density Up and tt-b. 
Then, Eqs. H25() and H26() can be cast into a simpler form. 
From Eq. (|26|l . we obtain the bosonic density tib = {^J■B — 
UpBnp)/UB and inserting this result into Eq. if^ . we 
obtain 



li[np\ = /ip 



Ub 



Ub 



(27) 
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Using H23() we arrive at an implicit equation for rip at 
fixed global chemical potentials /ip f-nd /ij,, 



"■F = J deN{e) f ^e-//p + -^/iB--^7 



(28) 



A solution of (|28|) with fermionic and bosonic densities 
rip and rie is a local minimum of the functional (|21|l . if 
the Hessian 



H 



9„ppi[np] [/pB 



(29) 



is positive definite. The derivative 9„p/i[np] is related to 
the response function x{T, q) via the compressibility sum 
rule and takes the form 



9„pAiK] = deNie)[-dJ{e-fi)]j (30) 

= ~[xiT,0)r'- (31) 

A positive definite Hessian imposes the stability condi- 
tions 

TrH = |x(T,0)| + f/B >0, (32) 
detH = |x(T, 0)1-1 Ub - Ul^ > 0. (33) 

The first condition is always satisfied for a repulsive in- 
teraction J7b > between the bosons. For a regular den- 
sity of states N{£), the second stability condition at zero 
temperature T — becomes 



C/b 



Niep) = 1 



(34) 



and coincides with the condition for phase separation in 
Refs.inilll ep denotes the Fermi energy. 

However, for the system considered here with 2D 
fermions on a square lattice the density of states ex- 
hibits a van Hove singularity. A setup with a fermionic 
density rip matching the van Hove singularity at eo and 
arbitrary bosonic density rig then is of special interest. 
These densities have to satisfy the conditions H25|l and 
(|26|l . from which we find the appropriate global chemical 
potentials /J^ = C/pb^If + Ugfi-B and, using /Lt[np] = eoi 
/Ip = [/pB«B + eo- Fixing the global chemical potentials 
/Ip and Jig rather than the densities and Ug, we can 
allow for the appearance of new solutions, i.e., phase sep- 
aration. The new fermion density rip then has to satisfy 
the condition (|28|l at fixed values Jig and /Ip of the chem- 
ical potentials, 



deNie) f 



C/2 



np) 



(35) 



The self-consistency equation (|35|) has the obvious so- 
lution Tip — Up. This solution turns unstable at the 
point where the second stability condition 1)33(1 is vio- 
lated: defining the deviation Sn^ ~ and expand- 
ing (|35() around rip wc find that 



a{T) 6rip+l3{T){6rip)^ +j{T) (,5np)^ 



0, 



(36) 



where a(r) = [1 + (U^g/Ug) J deNie + eo)dJie)] « 
ao(T-rps), while /3(T) and 7(r) > depend only weakly 
on temperature. For a system with particle-hole symme- 
try the quadratic term vanishes, i.e., P(T) = 0; this is 
the case for the situation discussed above with the dis- 
persion relation H12|) and a half-filled band with Fermi 
energy ep = 0. The sign change in a(T) at Tps matches 
up with the violation of H33|) and the expression for the 
transition temperature Tps agrees with the one derived 
previously in (^0)1 . At high temperatures T > Tps, a > 
and Eq. (|35|l has only the trivial solution np = 0. For 
low temperatures T < Tps, a < and two new solutions 
rip ± Srip appear with 



drip 



v/^(Tps-r)i/2, 



(37) 



where we assume /3 = and T close to Tps- 

At zero temperature and for ApB <C 1 we find the den- 
sity change 



Srip = ^ 



[/2 

dnN { eo + -^n 



(38) 



(5npApB In 



/ 16c2 JpC/p 
\UUSnp\ 



with C2 = exp(l). Solving for Srip provides us with the 
following shifts in the fermionic and bosonic densities 



(39) 



... UgJp ( 1 

rip - np = onp = ±16c2— 2— exp - — 

f^pB V ApE 

U.n._ Jf / 1 

rig-Ug = — — dnp = Tl6c2— — exp — 

Ug Upg \ ApB 

The relation between the T = density shift Srip and the 
critical temperature Tpg takes the form 



=C2/ci«2.40. 

Ug Tps ^' 



(40) 



Inserting the solutions into the free energy pro- 
vides us with the energy shift per unit cell at zero tem- 
perature 

^ = -EkiSnpf. (41) 



C. Phase diagram 

The solution discussed above allows us to find the 
phase diagram of the system. At low temperature T < 
Tps, we distinguish between a low-density 'gas' phase 
with np < np — i5np, and a high-density 'fluid' phase 
with np > np -|- i5np. The bosonic density derives from 
rig = [jig — UpgUp) /Ug. For a repulsive interaction 
Upg > the fermionic low-density phase partners up 
with a high bosonic density, while for attractive inter- 
action Upg < the fermionic low-density phase implies a 
low bosonic density. 
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5np(r) 




• |If>^f 



. I^F<I^F 



PS 



FIG. 1: Uf-T phase diagram. The grey region de- 
notes the phase separated regime, while the sohd Une de- 
notes the Une of fixed chemical potential Jl-p. The dotted 
lines describe fixed chemical potentials {fi-p — Jlp)/Tps = 
±0.02,0.06, O.Of, 0.14, 0.18. The lines derive from Eq. ^ 
with A = 0.2 and the dispersion relation 1121 . 



The phase transition is of the standard liquid-gas tran- 
sition type, see Fig. ^ At temperatures below critical, 
T < Tps, the low-density 'gas' phase is separated by a 
first-order phase transition from the high-density 'liquid' 
phase. The transition is driven by the chemical potential 
and takes place at the critical value 



(42) 



cf. Eq. I|27(l . where we have used the results eg = and 
?Tf = 1/2 valid for the dispersion 1)12(1 . A fixed averaged 
fcrmionic density between rTp — Suj, < < + Sn^^ 
is only realized via coexistence of the low density 'gas' 
phase and the high density 'liquid' phase. The first order 
transition terminates in a critical endpoint at the temper- 
ature Tps and density rip =np. Varying the temperature 
across the critical value Tpg along the isochore Up = Wp, 
we obtain a second-order phase transition between the 
homogeneous phase and the phase separated state, see 
Fig. ^ This transition appears for arbitrary weak cou- 
pling due to the enhanced fermionic density of states 
for 2D fermions. 



IV. FINITE TRAPPING POTENTIAL AND 
CONCLUSIONS 

We have shown that optical lattices strongly modify 
the behavior of Bose- Fermi mixtures via the appearance 
of van Hove singularities. The effect is most pronounced 
in two-dimensions, where the density of states diverges 
logarithmically, see Eq. Hll|l . Then, for a fermionic chem- 
ical potential matching the position of the van Hove 
singularity, i.e., fi = eo, the mixture of bosons and 
fermions undergoes a second order phase transition at 





FIG. 2: Sketch of the fermionic and bosonic density profiles 
in the weak coupling limit Afb <C 1. (a) The fermionic den- 
sity exhibits a jump in the density by 2Snp as the first or- 
der transition line is crossed. At the same position in space, 
also the bosonic density exhibits a jump. For attractive in- 
teraction UpB < between the bosons and the fermions the 
bosonic density increases, see (b), while for repulsive interac- 
tion UpB > the bosonic density decreases at the jump (c). 



the critical temperature Tps . In analogy to the standard 
'liquid-gas' gas transition, a fermionic density Up with 
Up — Srip < Up <np + Sup is unstable and inaccessible in 
an experiment for temperatures T < Tpg. The fermionic 
density profile in a harmonic trap can be derived from 
Eqs. (|25|l and H2t)|) . and its behavior at zero temperature 
and weak coupling is sketched in Fig. |21 The weak cou- 
pling ApB ^ 1 between the fermions and bosons makes 
sure that the density profile is only slightly modified com- 
pared to the noninteracting limit. As the local chemical 
potential fi{x) matches up with the van Hove singularity 
£0, the fcrmionic density exhibits a jump by 2Snp related 
to the crossing of the first-order transition line. This be- 
havior does not depend on the sign of the interaction UpB 
between the bosons and fermions. In turn, the bosonic 
density profile strongly differs for attractive and repulsive 
UpB- At the same position as the fermionic density profile 
exhibits a jump in the density, the bosonic density pro- 
file also exhibits a jump by — 2C/fb<^^f/C/b- The bosonic 
density increases similarly to the fermionic density pro- 
file for an attractive interaction UpB < between the 
bosons and fermions, while it is decreased for a repulsive 
interaction UpB > 0, see Fig.|5| The different behavior of 
bosons and fermions is a consequence of the asymmetric 
role played by the constituents: the transition is driven 
by the van Hove singularity in the fermionic density of 
states and is independent on the bosonic density. Note 
that the situation is quite different in the strong coupling 
limit ApB ~ 1. Then the van Hove singularity in the den- 
sity of states plays a minor role and the nature of the 
phase separated state is dominated by the character of 
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the interaction and the trapping potential , a situation 
previously discussed in Refs. [l^, [13; IH; Uji Hit 

We have shown that the presence of van Hove singular- 
ities in the fermionic density of states produces new and 
interesting features in the context of phase separation in 
a 2D Bose-Fermi mixture. At the same time, these van 
Hove singularities also enhance the instabilities driving 
other quantum phases competing with phase separation. 
E.g., the instability towards BCS-superconductivity ^3 



or the recent proposal for a supersolid phase are 
both driven by Fermi-surface nesting, usually involving 
a \n{eF/T) divergence which is enhanced to a [ln(eF/T)]^ 
singularity in the presence of a van Hove singularity. 2D 
Bose-Fermi mixtures then are promising candidates for 
the observation of such new quantum phases, provided 
they successfully compete against the tendency towards 
phase separation lurking at Tpg. 
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